Is spin-charge separation observable in a transport experiment? 
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We consider a one-dimensional chain consisting of an interacting area coupled to non-interacting 
leads. Within the area, interaction is mediated by a local on-site repulsion. Using real time evolution 
within the Density Matrix Renormalisation Group (DMRG) scheme, we study the dynamics of wave 
packets in this two-terminal transport setup. In contrast to previous works, where excitations were 
created by adding potentials to the Hamiltonian, we explicitly create left moving single particle 
excitations in the right lead as the starting condition. Our simulations show that such a transport 
setup allows for a clear detection of spin-charge separation using time-resolved spin-polarised density 
measurements. 
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Introduction 

Although spin-charge separation (SCS) is undisputed 
to happen in one-dimensional systems, its direct experi- 
mental observation and verification is a subject of many 
works over the past years. In one dimension, the Fermi 
liquid picture of non-interacting quasi-fermions is not ap- 
propriate Instead, in the paradigmatic Tomonaga- 
Luttinger model and in the prototypical Hubbard model 
collective bosonic spin and charge excitations are the low- 
energy excitations. These so-called spinous carry spin 
but no charge and so-called holons carry the charge but 
no spin and spinous and holons may propagate with in- 
dependent velocities. 

Up to now few experimental setups barely show di- 
rect evidence of separate spin and charge excitations ■ 
Beside tunnelling experiments e.g. momentum-conserved 
tunnelling on cleaved-edge overgrown GaAs/AlGaAs het- 
erostructures yl and single-point tunnelling into carbon 
nanotubes 0, there is angle-resolved photo emis- 
sion spectroscopy (ARPES) on quasi-one-dimensional 
SrCu02 material Q and on the organic conductor TTF- 
TCNQ 0, i which seem to show distinct spinon and 
holon branches in the spectral function. Numerical sim- 
ulations explaining the latter experiments were done by 
Benthien et. al. [0]. 

Here, we propose the following transport experiment: 
A one-dimensional interacting region is connected to non- 
interacting leads on both sides. In one of the leads a sin- 
gle electron with average momentum kg is injected into 
the system. The excitation, by passing through the inter- 
acting region, ends up in the other non-interacting lead, 
where a spin-resolved and time-resolved measurement of 
charge density is carried out. This setup poses the fol- 
lowing questions which we answer in this paper: If one 
injects an electron with definite momentum at some time 
Iq into an interacting system, where the separation of 
spin and charge is known to happen, what happens to 
the electronic excitation within the interacting region? 
And what kind of excitation emerges from the interact- 
ing area at a later time? Will we see distinct excitations 
in spin and charge density, will we find spin and charge 



density to be recombined to one full electronic excitation, 
or will we obtain an incoherent superposition of many ex- 
citations? 

Safi and Schulz [To| analyse such a transport setup ana- 
lytically for a spinless Tomonaga-Luttinger liquid model. 
For the spinful case, however, they have to neglect "umk- 
lapp processes and the backscattering of electrons with 
opposite spin" to sketch a qualitative picture, thus look- 
ing at an idealised system, where spin and charge currents 
are conserved separately. The first numerical simulations 
for spin-charge separation were performed by Jagla et. al. 
11 1 using exact diagonalisation for 16 sites. Zacher et. al. 
12| used Quantum Monte Carlo simulations to study the 
spin and charge susceptibilities of the Hubbard model. 
KoUath et. al. [l3| presented the first spin-charge sep- 
aration calculations in the framework of adaptive time 
dependent DMRG (td-DMRG). This calculation was re- 
done by Schmitteckert [l^l , keeping up to 5000 states per 
DMRG block showing that the adaptive time evolution 
scheme is prone to large numerical errors in the long time 
limit. In the same line, Schmitteckert and Schneider [Tsj 
used over 10000 states per DMRG block to show SCS in 
a 33 site system with periodic boundary conditions which 
avoids Friedel oscillations. 

In order to avoid the numerical problem associated 
with the adaptive time evolution scheme and to avoid 
the numerical costs of a pure full td-DMRG we apply a 
hybrid of both methods, see [l ^. [l6l . We first use the 



standard ground state DMRG [17| to obtain our initial 
state. We then perform typically ten time steps using 
the full td-DMRG jlSl and continue with an adaptive 
td-DMRG [ilJlO, |21| scheme. For recent reviews on 
DMRG, see [22 1 . In addition, two main objectives dis- 
tinguish our work from others that use real-time anal- 
ysis and model strongly interacting fermionic systems. 
1) Previous works either treated only homogeneous in- 
teracting systems [ill IHj 113 or modelled optical traps 
[l^. In contrast to that our system is a two-terminal 
setup modelling a transport experiment. 2) While in 
earlier works [ij, [3 excitations were studied based on 
the modification of local potentials, we create an excited 
state by explicitly adding one electron to the system and 
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we monitor the time evolution of such an excited state. 

Outline First we motivate our creation of an excita- 
tion and outhnc how the numerical method is used to 
measure the obscrvablcs. This is accompanied by bulk 
simulations on a non-interacting and an interacting sys- 
tem. Then we apply the method to the transport setup. 

Motivation and method 

An excitation, created at a given site xq in the tight- 
binding model by using cj.^, excites all eigenmodes in 
the diagonal basis of momentum eigenstates of the tight- 
binding model. This is not particularly useful since we 
want to consider electrons with a definite momentum. 
On the other hand, creating an excitation using c[.^ lo- 
cal in the momentum space is again an eigenstate of the 
system and invariant under time evolution. Therefore in 
the ground state of a half-filled tight-binding model 
(for each spin specie a = i^), 

Hth,a = - XI (4-i,<tCx,^ + h.c.) (1) 

X 

we create an excitation with spin a by using a Gaussian 
distribution of creation (annihilation) operators, with 
momenta centred around fco > kp^ko < kp): 

g(t)(fco) = Xe-'^ci^) (2) 
fc 

M 

x = l 

where k are the momentum eigenstates k € {— tt -|- 
. . . , tt}, and ctq is the width of the excitation in mo- 
mentum space. Additionally, the excitation will be cen- 
tred around Xq and normalised such that integrating over 
all k will resolve to exactly one added electron (hole) . A 
visualisation of this construction is shown in figure [TJ An 
excitation created in this way 

with C = (*| (5i^^(fco))'^5a^(fco) I*) will move through 
the system in one direction with a group velocity w(fco) = 
^Ofc I ko ■ Using the Schrodinger picture and evolving our 
initial state |\I'+^(t)) = e*^* j^'"''^), the measured observ- 
able is the spin-resolved density 

n,{x,t) = {^+\t)\cl,c,^,\^+\t)) (4) 

which allows the measurement of charge density ric = 
n| -|- ni and spin density Ug = ■^{ni — n^) at any time t 
as a function of x. 

In the following figures, we use DMRG [T^ to generate 
the ground state j'l') and the excited state |^^^)- For 




FIG. 1: (Colour online) Visualisation of the creation of an 
electron (left) or a hole (right) by indicating the occupation 
of single particle levels in the dispersion relation of the non- 
interacting problem. The dotted line represents the Fermi 
energy. 



both, the full and the adaptive td-DMRG, the time evolu- 
tion is performed by a Krylov subspace based evaluation 
of the time evolution \^+\t + At)) = e*^^* \^+\t)), 
for details see [H, [H, . The time steps were At = 0.5 
(in units of time) for both methods. 

We use hard-wall boundary conditions in all simula- 
tions and all energies arc in units of the hopping pa- 
rameter which is set to 1 in ([1]). Wc also calculate the 
static ground state density no^„{x) = cj. ,^c^ g. I^f) 
at time t = and form the corresponding charge and 
spin densities rto,c,"-o,s- Then we plot the quantities 
ns{x,t) — no^s{x) and nc{x,t) — riQ^cix). This trick evens 
out all stationary oscillations already present in the 
ground state, like the Friedel oscillations resulting from 
using hard- wall boundary conditions. Additional 2kF- 
oscillations resulting from the density distortion of the 
added excitation are evened out using a three-point av- 
erage. 

In our non-interacting (tight-binding) reference system 
^ in figure [5] we create an electron at the real space po- 
sition Xq = 15 (upper diagram) which travels to the right 
(lower diagram). The system of size M = 102 is at half- 
filling plus one electronic excitation N-\^ = 51-1-1. The 
Ni = 51 down-spins are, of course, irrelevant in the non- 
interacting case and constant in the time evolution. The 
excitation's momentum is centred around fcp = 7r/2 -|~ 0.1 
with width ctq = 0.03. With this width, we ensure that 
the momentum distribution is far away (compared to the 
width) from the Fermi surface fco ~ + 3(To but still as 
close as possible to the linear regime of the cosine band 
at kp- The snapshots at t = and t = 36 in figure 
[2] of spin and charge density serve as a proof-of-concept 
for the numerics. They show the expected synchronous 
motion of the electronic excitation with group velocity 
vg = 1.93 ± 0.02. A delta-function-shaped excitation 
in the thermodynamic limit of this system would have a 
group velocity of wg(^o) ^ 1.99. Since we are in a finite 
system, the discretisation and the cut-offs at the Fermi 
energy and the band edge directly influence the actual 
shape of the initial density distribution. Also, the non- 
linearity of the band affects the packet propagation via 
the width and position of the momentum distribution. 
The (otherwise perfect time invariant) Gaussian shape is 
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FIG. 2: (Colour online) Spin (solid, red) and charge (dashed, 
blue) densities of a td-DMRG simulation of one additional 
up-electron over the background of a half-filled tight-binding 
model for each spin specie. The Gaussian excitation on the 
left at time t = (upper diagram) moving to the right at t — 
36 (lower diagram) retains its shape. Averaging is described 
in the text. 



distorted and the packet is slowed down. This system was 
calculated using A^cut = 1000 states per DMRG block, 
applying 10 time steps in the full td-DMRG and keep- 
ing 1000 states for the consecutive adaptive td-DMRG. 
The discarded entropy never exceeds the order of 10"^. 
Comparing to an exact diagonalisation (ED) we define 
the peak relative error in the densities by 

\Sn,it)\ ^ max ( I^EhI^iII^™^!^) (5) 

This quantity reaches from initial |(5rip(0)| = 7-10~^ peak 
relative error to |5ri,p(36)| = 7 • 10""^ peak relative error 
a.t t — 36. 

Now we repeat the simulation in a Hubbard model bulk 
system 



bulk 



E 



(6) 

While the previous excitation was a superposition of well 
defined eigenstates in the tight-binding model, it is not 
clear which eigenmodes are excited by adding/removing 
an electron from an interacting system. Nevertheless, 
using the same cite in equation [3] to create an excitation 
still results in a Gaussian shape in spin and charge den- 
sities (see figure |3]) , although the corresponding Fourier 

transformed operators c^^^ '^o longer create single par- 
ticle eigenstates of the system. This time we insert an 
excitation corresponding to one (up-spin) hole into the 
ground state of a repulsive particle-hole symmetric Hub- 
bard model dH) with t/ = 4 at a fiUing of |i w 0.46 with 
a central momentum fco — 7r/2 — 0.5 ~ kp — 7(To and a 



width of ctq 0.053. The calculations were carried out 
with A'^cut = 2500 states kept per DMRG block in the 
full and adaptive time dependent DMRG. One clearly 
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FIG. 3: (Colour online) Spin and charge density (see legend 
of fig. (2)1 in a Hubbard model at time t — (upper diagram) 
and t = 25 (lower diagram) show spin-charge separation. This 
is data from a td-DMRG simulation on 90 sites system with 
A''| = 41 — 1, A'^l = 41 electrons and a Hubbard interaction of 
f/ = 4. 

observes SCS driven by the Hubbard interaction between 
t = (upper diagram) and t = 25 (lower diagram). We 
extract the group velocities from the main peaks to be 
vs = 1.12 ± 0.05 and vc 2.03 ± 0.05. The velocity is 
constant through out the simulation and the error stems 
from the read-off error. The one dimensional Hubbard 
model can be solved exactly |2^. From the analytic ex- 
pression by Coll [l^], the spin and charge velocity of the 
low energy spinous and holons can be derived as done, 
e.g. by Schulz (25j . The excitation's momentum of 
corresponds to a band filling of n = 0.68 in units of the 
hopping. Given the error estimate above, it is accurate 
enough to read off a charge velocity of vc.ono hoion ^ 2.05 
and Vs. one spinon ^1.15 (cf. figurc 11 of the first paper of 
Rcf. [2^). Although we are with 90 sites far from sim- 
ulating the thermodynamic limit, the group velocities of 
our spin and charge density fit remarkably well to the 
Bcthe Ansatz results for the velocity of one spinon and 
one holon. A similar good agreement was found for a 
spinless Luttinger liquid simulation [T§ |. 

Our electronic excitation with a definite momentum 
fco created a unidirectional running peak in the non- 
interacting system. Here, however, our initial excitation 
creates both a right and a left moving density peak in the 
charge sector. At time t = 25, the right moving charge 
peak has aheady been reflected on the right wall and 
is located at around x = 60. The spin density, in con- 
trast, displays predominantely a left moving peak. Fi- 
nally, there is charge density accumulation trailing the 
main charge peak and at the location of the spin peak. 
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These might be more compheated features of the interac- 
tion or result from the finiteness of the system as in the 
non-interacting example. 

Following the previous line of arguments, it seems rea- 
sonable to do a transport setup that uses non-interacting 
leads as a starting point. We let this excitation travel 
through an interacting region where SCS occurs and 
then let the resulting excitations travel again into a non- 
interacting area. In the starting lead our electronic ex- 
citation is well understood. When this excitation travels 
into the interacting region, we know the momentum and 
energy of the injected wave packet (one electron or one 
hole) but the questions arise: what kind of excitation 
leaves the structure and what can be measured on the 
opposite (non-interacting) lead? 
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FIG. 4: (Colour online) Excitation in a transport experiment: 
a created hole in the right lead (upper diagram, at t = 0) 
passes the interacting nanostructure (black bar) undergoing 
SCS (central diagram, at t = 17). At t = 35,40 (lower di- 
agram) spin and charge densities travel independently with 
equal velocity in the left lead. 

We use a system of M = 100 sites, split into 41, 29 and 
30 sites for the different areas and a total of 7V| =48—1, 
Ni = 48 electrons. The corresponding Hamiltonian reads 

-f^transport — ^ Hth,(j + V^ate ^ {^x.^ + nx,l) 
cr=T,J. 31=42. .70 

where i?tb,o- was defined in ([1]) and f/ = 2.5. A gate 
voltage on the interacting area levels the Fermi surface 



to a half-integer number of up and down electrons on 
the structure {N-^ = iV| = 12.5) and minimises reflec- 
tions into and out of the area (see figure |4|) . This lo- 
cal chemical potential V^atc in Q was pre-calculated in 
a self-consistent way. Thus we have a Hubbard model 
at ~ 0.43-filling in the central (black bar) area and 
half-filled tight-binding leads. We also chose parameters 
fco — 0.437r — 2(70, "'o = 0.03 to maximise tunnelling while 
staying as close as possible to a quasi-linear dispersion. 
The number of states kept was A*"cut = 2500. Figure [3] 
shows the time evolution of spin and charge densities for 
several times. There are reflections of charge density at 
both edges of the structure. They are suppressed by the 
choice of the chemical potential as discussed above. The 
majority of spin and charge density transmits through 
the nanostructure. SCS clearly happens inside the in- 
teracting region (central diagram at t ~ 17) and SCS 
remains intact after leaving the interacting area into the 
opposite lead (lower diagram at i = 35 and t = 40). Here 
spin and charge now travel along separately but again 
with identical constant velocity. 



Conclusions 

In our simulation of a microscopic experiment on an in- 
teracting nanostructure we saw that the spin and charge 
excitation of an injected electron can be separated as it 
is known to happen in one-dimensional interacting sys- 
tems. The spin and charge separation can be directly 
observed by a time-resolved measurement of the spin- 
polarised density. There is no need for spin and charge 
to recombine to a complete electron (hole) for a measure- 
ment of a single electron (hole). Note that the charge and 
spin excitations in the non-interacting leads are valid ex- 
citations of the Fermi liquid. They are not elementary 
excitations of the Fermi liquid, instead they are a com- 
position of (at least) two excitations of the Fermi liquid. 
As an example, think of the combinations {c\ ^ ± c\ 
However the resulting distribution of separated densities 
originated from a single particle injection into the Fermi 
sea. Thereby, we confirm for the Hubbard model, what 
Safi and Schulz [l^ sketched for a low-energy Luttinger 
liquid transport setup with restricted interaction. It is 
interesting that our small system already shows the fin- 
gerprint of Luttinger Liquid description expected in the 
thermodynamic limit. Also note that even for highly en- 
ergetic incoming electrons we find the spin-charge sep- 
aration, albeit obscured by band effects. With the td- 
DMRG method even more realistic models, like the ex- 
tended Hubbard model with a finite next-nearest neigh- 
bour interaction could be investigated. 

Experimentally, the measurement of the densities 
would have to happen on a time scale which is set by 
the Fermi velocity and the length of the interacting sys- 
tem. For nanoscaled systems and metallic excitation ve- 
locities tens of Terahertz resolution would be required, 
which might be achieved using optical techniques. Low- 
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ering the Fermi, holon and spinon velocities would en- 
hance the detection possibilities. It is of general interest 
to know, whether a series of quantum dots or any other 
strongly correlated system with a much smaller Fermi 
velocity can provide the ingredients for the realisation of 
such a transport setup. 
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